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ABSTRACT 


Inertial waves excited in the mixed layer by hurricane 
Frederic, had hcrizontal scales of approximately 1 to 2 
times *he tarcclinic Rossty radius of deformation (50 km) of 
the first mode near the Lesoto Canyon. Initially, energy 
propagated vertically at akout 1.25 km/d and horizontally at 
80 kmyd. These waves spun down over e-folding scaies of 
four inertial periods as energy propagated vertically at 270 
m/d and horizcntally at 30 kn/d. Inertio-gravity waves in 
the deep thersccline had hcrizontal scales of 25 to 50 km 
and vertical scales approximately equal to the water depth. 
The energy of these waves was dcminated by the barotropic 
mode with scme contributicns frem modes 1 and 2. These 
waves were not admitted to the shelf region because the bct- 
tom Slope was greater than the slcpe of the internal wave 
characteristics. 

The mean flew followed the  iscbaths at all levels, tut 
it was in the cpposite direction in the bottom layer. The 
mean flow initially decreased alcng the eastern boundary of 
the canyon as the storm fercing readjusted the flow. Near- 
bottcm temperature variations of 49C were associated with 
the storm surge and advection in the along-track direction, 


particularly along the north rim of the canyon. 
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I. INTRODUCTION 


" +he doer alone learneth " 


Nietzsche >n. 


The oceanic response tc intense, transient atmospheric 
events, such as hurricanes, is dcminated by robust inertial 
wave excitaticn in the mixed layer and the subsequent propa- 
ga+icn and dispersion of inertio-gravity waves in the ther- 
mocline. The study of these forced waves has been largely 
*through numerical and analytical investigations, because of 
the lack of sufficient data to resolve the scales and ener- 
getics of the motion. Hence, observations of ocean currents 
and temperatures during a hurricane are neccessary for vali- 
dating models and thecries dealing with forced inertio-grav- 
ity waves. 

A comprehensive data set of ocean observations vas col- 
lected by the U. S. Naval Oceanographic Office (NAVOCEANO) 
during the passage of hurricane Frederic in the summer of 
1979 (Shay and Tanul, 1980). NAVOCEANO deployed three 
moored taut-wire current meter arrays in the northern Gulf 
of Mexico soutkeast of Mobile, Alabama. During the period 


of deployment, hurricane Frederic passed within 80 to 130 kn 
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77 "he array sites (Fig. 1) at 2100 GMT 12 September. 


Five 
hcurs later, Frederic made landfall at Dauphine Island, Ala- 
bama. During and subsequent to the passage of Frederic, the 


Naticnal Hurricane Research Division dropped several expend- 


able bathythermcgraphs (AXET) from reconnassiance aircraft 


j-- 


n the area cf the current meter deployments (Black, 1983). 
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Figure 1. Back of Ene Frederic. The contours are in 
thems, 5570 5 iticns of current meter 

rays, Ad EANO) depicts the position of 

ET drops (Black, 1983). 

e 


The track of hurricane 
deric is based on a 000 by Hebert (1979). 
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This thesis addresses the problen of the vertical and 
horizental dispersion of the fcrced inertial wave energy 
from the mixed layer. The forcing readjusts the mean flow 
as well as generating anisctropic inertial motion. The cur- 
rent meter observations allow determination of the scales 
and the energetics of the forced, anisotropic wave motion 
<hrough use cf spectrum analysis and complex demodulation. 
The modal structure is examined by formulating the Sturn- 
Liouville problem and then performing a least squares fit 
between the eigenfunctiors and the demodulated time series 


Sethe current observations. 
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II. HISTORICAL REVIEW 


A nonlinear model for a stationary hurricane (O'Brien 
and Reid, 1967) was developed to simulate the observations 
made in the wake of hurricane Hilda (Leipper, 1967). The 
results indicated that upwelling regions were restricted to 
the area enclcsed by the maximum wind regime, whereas down- 
welling occurred outside this regime as warm water was 
advected away from the stcrm center. 

Cne of tke most proncunced responses of the ocean to 
wind forcing is the generation of inertial oscillations in 
the surface mixed layer. Polla rd (1970) simulated the gen- 
eraticn of these waves using a two-dimensional, vind-driven 
model developed by Pcllard and Millard (1970). The pre- 
dicted amrplituces and decay rates cf the forced oscillations 
agreed well with the Woods Hole Oceanographic Institution 
(WHOI) 'site D' cbservaticns. Using the same data set, Pol- 
lard (1980) showed that under relatively strong wind condi- 
tions, 67 % of the horizcntal kinetic energy (HKE) in the 
mixed layer was found near the inertial frequency. The sub- 


sequent dcwnward propagation of energy was estimated to be 
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7” “he order cf 10-3 cm/s ( 1 nd). Ties vertical group 
velocity was too small to account for the radiation of iner- 
+ial waves and the subsequent loss of energy from the mixed 
layer. The cerresponding scales cf the inertial wave motion 
were about 100 “o 240 m in the vertical and hundreds of 


kilometers in the horizontal (Table I). 


TABLE I 
Compariscn of Inertial-Internal Wave Parameters 


Frcm Brocks (1983) 


| Length Brooks Price Pollard | 
| (1983) (1983) (1980) | 
| Scales | 
Hcrizon*al(km) 370 480 700-1700 
| Vertical(n) 1000 1000 100-240 | 
| 
Group. | 
| city | 
| Herizontal (km 7d) 23 86 0.8-17 | 
| Vertical (1/d) 60 160 0.03-3 | 


Using a linear, two- layer ncdel, Geisler (1970) simu- 
lated inertic-gravity waves as part of the baroclinic 
response in the wake of a translating storm. Tbe critical 
factcrs governing the generation of these waves were: 


e the translational speed of the hurricane must be 
greater tkan the internal wave phase speed; and 


ә the hcrizcntal scales of these waves must be comparable 
to the oceanic Ressby radius of deformation. 
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He furzher noted that the ocean response to a moving hurri- 
cane is strcngly baroclinic. 

The ocean thermal response to hurricane forcing, was 
investigated numerically ty imposing Ekman layer dynamics in 
a mixed layer model (Elsberry et al., 1976). Фиситари от 
energy frcm the wind stress was forced to gen¢rate inertial 
oscillations and cause turbulent mixing via entrainment. 
Upwelling alsc enhances tke turbulent mixing process in the 
upper layers of the ocean, The main result was that advec- 
tion dominated tke thermal response near the storm track as 
opposed to the heat loss from the ocean surface to the 
storm. 

Strong atmospheric fcrcing enhances the turbulent mixinq 
process in the surface mixed layer and causes it to deepen. 
As the mixed layer continues to deepen, the thermocline 
begins to erode as water from the thermocline is entrained 
into the mixed layer. However, stratification tends to sup- 
press turbulent mixing and aids in the creation of small- 
scale internal waves either through entrainment or 
Kelvin-Helmholtz instabilites (KB) at the base of the mixed 


layer (Pollar et al., 1973 ; Garwocd, 1977). 
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Eelow the mixed Waye: ocean current variability is 
linked to vertically propagating wave groups of the internal 
wave field (Kase and Clbers, 1979). The actual mechanisms 
for the vertical transport of energy from the wind forced 
mixed layer tc the thermccline for the generation of large 
scale inertial-internal waves are not well understood. Per- 
haps the most effective way of creating these large scale 
inertial-internal waves is through Ekman suction. Krauss 
972 a,b, 1976) showed that a horizontally varying wind 
stress causes a mass transport 90 deg. to the right of the 
stress in the Ekman surface layer. This surface layer div- 
ergence is acccmpanied by an upward displacement of the iso- 
pycnals or upwelling of cccler water from below which tends 
to suppress the turbulent mixing process. During the relax- 
ation of the wind, these isopycnals are displaced downward , 
which contributes to the creaticn cf large scala inertial- 
internal waves at the base of the mixed layer. 

During the ¡passage of hurricane Belle (Mayer et al., 
1501), ocean current and temperature measurements vere 
acquired on the continental shelf of the Middle Atlantic 
Bight. Their analyses indicated that most of the inertial- 


internal wave energy was ccntained in the first mode at the 
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dee per sites (water depth of 70 m) and in a heavily damped 
second mode at the shallow sites (water depth of 50 n). The 
variability in the modes was not attributed to the spatial 
HBezlability of ee wind stress, The bottom slope ard mean 
velocity fields significantly altered the oceanic response 
to hurricane Eelle. 

A mixed layer model (Garwood, 1977) was embedded into a 
Multi-layer, primitive equation, ocean circulation model 
(GCM) by Adamec et al., (1980). Using this GCM model, Hop- 
kins (1982) simulated the baroclinic response to a forcing 
pattern similar to hurricane Frederic and compared the 
results to the data collected in the wake of Frederic (Shay 
and Tamul, 1980). The model predicted the inertial response 
in the mixed layer quite well. However, the predicted ocean 
current response in the sutsurface layers was less energetic 
than the ckbservations, and decayed much too fast. These 
discrepancies between the model simulations and observations 
meen due in part to the impesition of the rigid lid conditior 
at the surface, which eligzinates the barotropic mode. EH 
thermere, these model simulations did not include topogra- 
phy, even thcugh the observations are from a region where 
the kottom topcgraphy is rugged and the ocean depth is less 


than 1 kn deep. 
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Recently, ocean measurements vere made during the 
passage cf hurricane Allen in the western Gulf of Mexico 
rocks, 1983). The spatial scales of the inertial wave 
were different fron the 'site D' cbservations, as shown in 
Table I. The vertical scales in the ocean are generally 
much greater under hurricane forcing than during the passage 
27 cold front, while the converse is true for the horizon- 
tal scales cf motion. Dee verccalugsoup"veloc2i-y, which 
depends on Ecth the wave and the Brunt-Vaisala frequencies, 
for the Allen cbservations exceeded the group velocity from 
the "site D! observations. This estimate from the Allen 
observations is misleading due tc the lack of upper thermo- 
cline and mixed layer data. Moreover, the wind speed is 
much greater during the passage cf a hurricane than a fron- 
tal passage. Therefore, the amcunt of turbulen- mixing in 
the mixed layer and inertial wave excitation should be much 
more rapid during a hurricane. Other fea-ures measured in 
the Allen response included topcgraphical dependence, and 
the vertical phase locking during the first few inertial 
pericd (IP) fcllcwing the storn. 

The baroclinic response of the ocean to a hurricane was 


modeled by (Erice, 1983) using a multi-level, inviscid 
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model. The scales and energetics of the simulated inertial 
respcnse were similar to mixed layer data collected by a 
NOAA data kuoy during the passage of hurricane Eloise. The 
horizontal and vertical scales, as given in Table I, vwere 
quite large in ccmparison to the thickness of the thermo- 
cline (200 n). The maxirun energy of the inertial-internal 
wave moticn vas predicted to occur at a distance of twice 
the radius of maximum winds, hereafter referred to as the 
maximum wind regime, which was roughly 80 km for Eloise. 
Further, the rate of vertical energy propagation was large 
in ccmparison to both the Allen and ‘site D' observations, 
and accounted fcr the depletion of energy from the mixed 
layer. Other mechanisms, such as KH instability (Pollard e+ 
2: 1973 : Garwood, 1977) and turbulence (Bell, 1978), act 
to remove energy from the mixed layer. However, these phe- 
nomena were neglected as mixing was not explicitly included 
me che model forrulation. 

Greatbatch (1983) modeled the nonlinear response of the 
ocean to a mcving stcrn. The major results were that the 
transition between upwelling and downwelling zones in the 
oscillatory wake was rapid, and that nonlinearities account 
for the displacement of the maximum response to the right of 


the storm track. 
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III. DATA AND ANALYTICAL 


— = cx eme До — == 


A. CURRENT METER DATA 

Ten Aanderaa RCM-5 current meters were deployed on three 
moored taut-wire arrays ir depths ranging from 100 to 470 a 
in the northern Gulf of Mexico. These current meters san- 
pled ccean current speed, direction and temperatures at 10 
minute intervals. Two current meter arrays (CMA2, and 3) 
were deployed on adjacent sides of the DeSoto Canyon where 
bathymetric ccntcurs converge to form the head of the canyon 
(see Fig. 1). Ihe other ıcoring was deployed closer to the 
coast in a depth of about 100 m of water where the isobaths 
are nearly parallel to tke ccast. A Synopsis is given in 
Table II cf tte storm period observazions, which extends 
from a few days prior to hurricane passage (12 September) to 
the end of the deployment feriod (mid-October). 

The quality cf the data is generally good; however, the 
lengths of the time series are not all the same. For 
instance, tbe Savonius rctors were eventually lost from all 
current meters in the mixed layer due “0 the large current 


speeds. These large mixed layer current speeds are not 
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Berrec+ted fcr rctor pumping, Wiser es a function- of ths 
mooring design as well as the large direction vane of th2 
Aanderaa current metərs (Fofonoff and Ercan, 1967) . Tn 
addition tc tke rotor problems, temperature records from all 
mixed layer current meters were unrecoverable because the 
Ocean temperature exceeded the threshold temperature of the 
Bhermistors of 21.5 °C . Sea surface temperatures 325 km 
south-southeast of the mocring sites, as measured by NCAA 


Buoy 42003, reacted 28.89C (Johnson and Renwick, 1981). 


Be STORM TRACK 

The stcrm track is based on the best position data on 
hurricane Frederic's movement through the Gulf of Mexico 
Webert,197/9). Because of the intensity of Frederic, recon- 
naissance aircraft corstartly monitored the storm. Frederic 
increased to maximum strength, maximum winds or minimum 
barometric pressure, 80 to 130 kn west of the CMA sitss 
about 2100 GMT 12 Septenber (see Fig. 1). A visible photo- 
graph of Frederic from a GOES satellite at 2001 GMT 12 
September clearly delineates a well developed eye of about 
40 tc 50 km in diameter (Figure 2). The translational speed 
of the hurricane at this time was about 7 to 7.5 m/s as it 


apprcached the Gulf Coast, which was much larger than the 
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internal wave phase speed. Thus, an inertio-gravity wavs 
response is expected in the wake of hurricane Frederic 


(Geisler, 1970). 


TABLE II 
A Syncpsis of Stcrm Period Observations 


Meter Reccrd Start End 
Depth Length Time Time 
(a) (days) (GMT) (GMT) Variables 
CHAT 
21 0 Nc hurricane data 
49 22:21 0040 5 Sep. 0600 27 Sep. V 
64 16. 76 0010 5 Sep. 1920 21 Sep. va E 
H2. X 252557 0040 5 Sep. 2040 25 Sep. n 
CMAZ 
19 24. €7 0040 5 1650 29 Sep. u,v 
179 2 1920 7 Sep. 0100 30 Sep. unu. 
324 24.76 0040 5 ə 1900 29 Sep. Wey T 
CMA3 
21 24.91 0100 2 Sep 2200 26 Sep. u,v 
251 36.00 0100 2 Sep. 0100 38 Oct. PAEA 
437 35:69 0100 2 Sep. 1730 7 Oct. AT T 
457 21.4 01C0 2 Sep. 9000 4 Oct. uz. T 
* time clcck synchronization problems 
u,v = horizontal velocity components 
T = temperature i 
GMT = Greenwich Mean Tine 
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Figure 2. Satellite Photograph cf Ereder 6 700 1:6451:12 
Sept. 1979. Visual satellite imagery is courtesy 


of NOAA/NESDIS. 
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5 WIND FIELD 

Wind field data were obtained from shore stations at 
Mobile, Alabama and Pensacola, Florida. These records indi- 
cated that the surface wind speed did not exceed 40 mys. 
However, cther reports suggested that wind speeds ranged 
between 48 +0 58 m/s as Frederic made landfall (Hebert, 
1979). These discrepancies are attributed to local boundary 
effects and the non-representativeness of wind data  ccl- 
lected in the coastal regicn. Marine winds vere also meas- 
ured Гу NOAA Euoy 42003 (Jchnson and Renwick, 1981). Wind 
speeds at this location never exceeded 35 m/s. The eye 
clearly passed over the bucy as indicated by the minimun in 
the wind field as the direction changed from 40 *o 200 9 
True (Fig. 3), with a corresponding decrease in prsssure to 


959 millibars (mb). 


D. VERTICAL TEMEERATURE GRADIENTS 

The AXBT data collected by the National Hurricane 
Research Divisicn, as reported by Black (1983), are used for 
the ccmputaticn of the Brunt-Vaisala frequency. These data 
were collected after the passage of Frederic in the area of 
current meter array deployments near the DeSoto Canyon (see 


7. 1. 


27 





36 


32 


3 
О 28 
= 
Ti 24 
ui 
o. 
Y) 
a 
z 20 
= 
18 
12 
8 
12 15 
ui 
ЫЈ 240 
ac 
= 
(2 200 
ЫЈ 
с 
(D 
w 160 
a 
= 
О 120 
- 
Q 
ЫЈ 
C 20 
a 
a 
È 40 
3 
0 
12 15 
Figure 3. 


-— 
00 
N 
=ә 
о 
ә 
О 
о 
-— 
M 
—_ 
m 
= 
00 
N 
c 
o 
Q 
o) 
o 
=- 
N 
— 
m 


18 21 0 3 6 9 12 15 18 21 0 3 6 9 12 15 
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JULIAN DAY/TIME (GMT) 
Wind Speed and Direction Measured by NOAA Buoy 


42003. The upper panel represents the observed 
wind speed time series where the abscissa depicts 
time in Julian Days starting on 1200 GMT 11 Sept. 
to 1500 GMT 13 Sept. 1979. Tha lower panel 
represents the concurrent wind direction time 
series. The hatched area depicts the period when 
the data buoy was in the eye of the hurricane 
(Johnson and Renwick, 19817. 





The vertical density profiles are computed at 10 m 
intervals using the AXBT data and climatological data from 


NAVOCEANO thrcugh the equation of state 


р = p,(l-aT-T)) 7 (1) 


where p is the density, T is the cbserved temperature, To 
and p. are tke reference temperature and density respec- 
mively frcm climatology, and a is the thermal expansion 
coefficient, which is taken to be 0.0002/9C . Some error is 
expected due tc the neglect of the salinity term. However, 
BwSwamina*ion of the clisatological T-S diagram, for the 
DeSoto Canyon area indicates that density variations are due 
to temperature rather thar salinity effects. The Brunt-Vai- 
sala, N? , frequency is computed using a certered finitə 


difference frcm the expression 


iz ? (2) 


2 
| 
! 
D рә 
D» 


O 


where g is the acceleraticr due to gravity andAp/Az is the 
Nemtical density gradient. Since the AXBT “data only 
extended to about 250 m, it was necessary to extrapolate the 


vertical temperature gradient to 470 m. The temperature is 
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assumed to decrease uniformly at a rate of 0.0311 °C/m from 
250 to 470 m (Fig. 4). -: The corresponding uniform increase 


mor the density produces a constant Erunt-Vaisala profile. 
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E. TEMPORAL AND SPATIAL SCALES 
The lccal variations in the bottom topography at the 


DeSotc Canyon (see Fig. 1) are appproximately parallel and 
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normal to the coast at CMAz and CMA3. Por example, Torth- 


ck 


Euuth flow a* CMA2 is in the cress-shelf direction, but it 
Hon the aleng-shelf direction at CMA3. Henceforth, a 
Cartesian coordinate system is used to represent the hori- 
zontal velccity components in the x and y directions, where 
moticn to the east and tte north corresponds to +u and +v 
current velccities, respectively. 

The fundamental time scale used in this thesis is an 
inertial ретјсб, which is equal to 24.10, 24.25 and 24.47 h 
at CMA 1,2 and 2 respectively. mehe following calcula- 
tions, the Inertial Period (IP) will be set equal “о 24.4 h. 
The intrinsic length scales are either the radius of maximum 
winds, approximately 40 kg for Frederic, or the Rossby radii 
of deformaticn, which are given ir Table III for the baro- 
tropic and first three barcclinic modes. These deformation 
radii are kased cn the inertio-gravity wavs phase speed com- 
puted from the Sturm-Liouville problen, which will be dis- 
cussed in a later section. For convenience, the scales in 
the x and y directions will Һе referenced as cross and 
along-track scales respectively, because the storm track was 


nearly perpendicular to the coastline. 


31 





TABLE III 


Rossby Radii of Deformation 


Phase Deformation 
Mode Speed Radili 
No. (m/s) (ka) 
0 68.0 950. 
1 4.0 56. 
2 2. 2 30. 
3 1.3 18. 


p EDMPUTATION CF SPECTRA 

The energy and cross spectra are computed in the follow- 
ing manner: (1) the mean is computed and is extracted from 
the time series; (2) a Tukey data window is applied to the 
data; (3) the data are transformed into frequency space via 
the Fast Fourier Transform (FFT); and (4) the transformed 
data are spectrally averaged over bandwidths which vary in 
Width with frequency, andare closely tied to the computa- 
tion of confidence levels. Тһе number of data points trans- 


formed by the FFT used in the spectrun analyses does not 


Ф 


have to be proportional to 20 , where n is aa integer. dh 
general FFT is advantageous in the analysis of events 
because of their transient nature. However, in the analysis 
of the spectra, equal numbers of data points are used to 


avoid leakage and smearing of energy fron adjacent frequency 
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bands (Otnes and Enochson, 1978). Finally, the higher fre- 
quency energy abcve the Nyquist frequency, 1/(2 A ¢) where 
A t is the sampling interval, is eliminated to prevent ali- 
asing of the spectra. 

The HKE is decomposed into the clockwis2 (CW) and coun- 
terclcckwise (CCW) rotating ccmponents using the rotary 
spectrum methcd outlined in Gonella (1972) and Mooers 
897 3). The main results of this analysis are “hə rotary 
spectrum, ellipse stability and orientation, and rotary 
coefficient. The orientation of the allipse indicates the 
direction cf the horizontal phase velocity, while the sta- 
bility indicates the variability in the ellipse orientation. 
A high stability index indicates that the wave motion is 
£ octropic, cz that the herizontal direction of phase pro- 
pagation is unidirectional. Conversely, a low index of sta- 
bility implies that the directicn of phase propagation is 
zandcm or isctropic. The Trotary coefficiənt indicates the 
amount of confidence that can be placed in the rotary spec- 
tral estimate and the rotational component that dominates 
the rotary spectrum. For instance, a +1.00 implies that the 
messen is polarized in the CW direction. A more rigorous 
mathematical treatment is given by Gonella (1972) and Mooers 


uo 7 2). 
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G. COMPLEX DEFODULATION 

A useful method for determin ing inertial oscillations in 
a time series of ocean ctrrent measurements is through the 
use cf the ccmplex de modulation method (Perkins, 1970). The 
mathematical details are given in Appendix A. This method 
can ke implesented either ky linear filtering or by perform- 
ing a least a e arno analysis on the time series. 
The linear filtering method is preferred because data 
points are lost whenever karmonic analyses are used to form 
the time-varying amplitudes and phases of frequency depen- 
dent moticns. This is especially useful in the analysis of 
events when all data points are needed to resolve che tran- 
sient moticns. A major disadvantage of the linear filtering 
2 cd is that it only applies to discrete frequencies, 
whereas the least squares method can treat all frequencies 
within the constraints of the time series. 

Initially, the data are band-pass filtered between 18 
and 30h to limit the frequency content in the time series 
to near-inertíal motions. These filtered data are multi- 
plied by the trigonometric arguments of the inertial fre- 
quency. At this point, some higher frequency motions hava 


been introduced into the time series (Appendix A) and these 
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are removed fy low-pass filtering the data at 22 h. The 
resultant time series is manipulated to form the instantane- 
ous amplitudes and phases cf the inertial motion. The final 
set of time-varying amplitudes and phases provide a local 
hazmcnic analysis, Gr more explicitly, a time series of 
spectral estimates for inertial motion. The filter used in 
the demodulaticn was a Lanzcos square taper window with 241 
weights (Fig. 5). 

The band-pass filtered, mixed layer velocity components 
at CMA3 are excellent examples for usin conplex demcdula- 
Son (Fig. 6). As the current velocities increase in 
respcnse to the storm, the amplitude of the forced inertial 
wave increases as well. Ihe iner+ial wave amplitudes modu- 
late the increases and decreases in the mixed layer velcci- 
Eros. After the maximum velocities occur, roughly 2 IP 
follcwing hurricane passage, the strength of the currents 
decreases linearly with e-folding scales of about 4 IP. 
This type of ¢-folding behavior is very commen in amplitude- 


modulated signals. 
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Figure 6.  Band-Fass Filtered Mixed layer V Velocity 
Ccmponent at CMA3. 


36 





17: TORM PERIOD 


In the fcllowing analyses, the data are divided into two 
pericds: 1) tte storm period observations as previously 
described in Table II; and 2) the quiescent period which 
encompasses the period frem the start of the data to a few 
days prior tc the storm. This chapter focuses on the iner- 
tial, subinertial and superinertial response in terms of 
energetics and scales of section to hurricane Frederic. The 
quiescent pericd observations and results are given in 
Appendix B. For convenience, the terns shallow and deep 
will refer tc the regicns where CMA1 and CMA2,3 were 


de plc yed. 


A. IİNERTIAL RESFCNSE 

Burricane Frederic passed abcut 80 kn vest of the DeSoto 
Canycn at 2100 GMT 12 September (UD 255) 1979. тһе ит те 
speeds began to increase at 0700 GMT or approximately 14 h 
PETOT to the time of closest approach. The currents in the 
mixed layer at CMA2 increased +o a maximum within an IP fol- 
lowing the passage of Frederic (Fig. 7). The storm caused 
the horizental velocity to oscillate with periods of about 


meron and to rotate in a CW direction, which indicates the 
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presence of inertial waves in the mixed layer. These wavss 
then decreased in amplitude over e-folding scales of 4 IP. 
The storm alsc caused the u-component to flow in the oppo- 
site direction from that cf the pre-storm flcw, whereas the 
v-component increased markedly in amplitude. The near-bct- 
tom current sreeds (shown later) began to increase within 3 


h after the mixed layer response. 
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Figure 7. Mixed layer Velccity Components at CMA2. The 
crdinate depicts the magnitude of the horizontal 
velccity components in cm/s for the east-west 
current (upper) and north-south current (lower) 
The hatched area depicts the period of hurrican 
passage. 
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In the mixed layer at CMA3, the storm amplified the 
pre-storm currents to a maximum of 135 cm/s (Fig. 8). AS at 
CMA2, this increase in the currents is a manifestation of 
inertial wave excitation. However, “he inertial response at 
CMA3 was different from CMA2 in that the kinetic energy of 
the waves was larger and persisted longer (approximately 2 
IP). The energy of these waves then decreased over e-fold- 
ing scales of 4 IP. The increase in near-bottom current 
(shown later) cccurred within 4.5 h following the mixed 
layer response. 

These hcrizontal differences of the inertial response in 
the mixed layer between CMA2 and CMA3 are observed in the 
amplitude changes of the HKE (Fig. 9). The initial increase 
in kinetic energy at the twc arrays was Similac. The naxi- 
mum in kinetic energy occurred first at CMA2 and led the 
mixed layer maximum at CMA3 by about 3 h. The differences 
in HKE levels were quite proncunced. For example, the mixed 
layer HKE level was roughly 60 % cf that found in the mixed 
layer at CMA3, even though CMA2 was closer to the eye of the 
storm. This large difference over a distance of 34 km may 
be due in part tc the strcnger mean current at CMA2, which 


may have reduced the amplitude of the inertial waves. The 
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Figure 8. Sisilar to Fig. 7, except at CMA3. 


direction of tke mean flow and the horizontal phase propaga- 
tion cf the inertial waves (shown below) roughly coincided, 
which sugges-s that there may have been some interaction 
between then. Fcllowing these maxima, the energy decayed 
with an e-folding scale cf 4 IP at both arrays. Forced 
inertial waves in the mixed layer that vary in space and 
time are one mechanism to drive inertio-gravity waves in the 


subsurface layers (Krauss, 1972a,b). 
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The vertical differences in the inertial response were 
even more dramatic. In the deer thermocline (251 m) at 
CMA3, the increase in the currents was sensed within 2 h 
follcwing the mixed layer response (Fig. 10). The initial 
maximum cccurred within an IP (JD 256) after the winds 
relax, which was followed by inertial oscillations in the 
horizcntal currents. In the presence of stratification, the 
particle  mcticns are elliptical rather than circular as 
observed in tte mixed layer. The anplitude of inertio-grav- 
ity waves decreased slightly following the initial response. 
About 10 IP later, the current increased to a secondary max- 
imum that nearly exceeded the initial maximum. The u-velo- 
city essentially defined an amplitude-modulation envelope 
that collapsed and expanded in time and persisted for 21 IP 
follcwing the initial maxiıun. The near-bottom currents at 
CMA3 (shown later) behaved in a similar manner 2xcept that 
the secondary maximum occurs 6 IP following the storm. 

The vertical changes in the currents at CMA3 are even 
more apparent in the prcgressive vector diagrams (PVD), 
which are Lagrargian representations of the Eulerian meas- 
urements. That is, the ccntinucus current measurements are 


integrated with respect tc time to approximate the 
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Mam ectory of the water parcels at that point. The horizcn- 
tal velocity components in the mixed layer at CMA3 defi- 
nitely showed the classical, CW rotating vector associated 
with inertial waves (Fig. 11a). Prior to the storm, small 
scale, freely prcpagating inertial waves were present, and 
they were subsequently azplified in response to the storm. 
These forced waves were superposed on a southward mean flow 
that increased in magnitude in response to the forcing. The 


horizontal displacement over successive time intervals of 2 
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IP increased followirg stcrm passage. After the maximum 
current velocities near JD 258, the amplitude of these 


forced, neéar-circular inertial oscillations decreased. 
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Figure 11. Progressive Vector Diagrams at CMA3. The * 
r Rs a time span of 48 h 952 on JD 
245 for_a) mixed layer (21 nm), and b) the deep 
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merhestchermeeline (Eig. 11b), the storm amplified the 
freely prcpagating inertic-gravity waves. However, the 
dimensions cf tte ellipses are distorted because of the 
increases in the mean flcw between JD 257 and 261. This 
pericd was immediately fcllowed by decreases in the rean 
flow until JL 265, when the inertio-gravity response was 
more pronounced. The pericd of increasing and then decreas- 
ing mean flow ccrresponds to the approximate modulation 
envelcpe and lag period cf 10 IP between JD 256 and 266 
(Fig. 10). After JD 266, the modulation envelope started to 
collapse as tke inertio-gravity waves decayed in amplitude 
With an e-fclding scale of about 4 IP. 

The vertical propagaticn of the energy, which is a con- 
sequence cf the generaticn of inertio-gravity waves in «he 
thermccline, is indicated by the vertical differences in the 
HKE between the mixed layer and the thermocline (Fig. 12). 
In the mixed layer, the maximum HKE occurred approximately 
one IP following the storm and persisted over the nex* 2 IP. 
Subsequently, a relative maximum was reached in the deep 
thermocline akcut 2 h after the surface maximum, and it was 
only 5 to 8 % cf the maximum HKE in the mixed layer. Based 


only cn these initial phase differences between the velocity 


45 





conpcnents at the two levels, the vertical propagation of 
energy would te about 3 cm/s (2.6 km/d). The thermocline 
HKE then remained nearly constant over the next 4 *o 6 IP. 
On JD 262, tte thermocline HKE gradually began to increase 
as the vertically propagating energy fron the mixed layer 
accumulated. Using this 10 IP lag period, the vertical pro- 
pagation of energy is about 0.03 cm/s (26 m/d). The "rire 
stimate cf 2.6 km/d is about 1 tc 2 orders of magnitude 
greater than the Eloise and Allen observations (see Table 
B). However, “he second estimate of 26 m/d is much more 
comparable to the Allen estimates. The propagation vel 
energy will be addressed below using the dispersion relation 
of inertio-gravity waves. The phase propagation was clearly 
upward over the first 10 IF following the storm. Over the 
subsequent 10 IP, the proragation was not as clear as phases 
in the mixed layer became erratic, presumably due to the 
decreasing energy levels cr spectral broadening. 

The inertial peak dcainated the HKE spectra at both 
locations and at all depths. The change in the inertial HKE 
between the surface and bcttom layers was about an order of 
Magnitud2, while the energies between mid-depth and interme- 


diate levels did not change significantly at CMA3 (Fig. 13). 
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Althcugh the inertial peak was prominent, the variability in 
spectral shapes was due tc the differences in the background 
continuum cf internal waves which varies with depth (Fu, 
1980) and Eandwidth averaging. 

The HKE srectra are decomposed into CW and CCW rotating 
compcnents to clarify the rotational characteristics of the 
inertial wave motion. Tke spectral peak near the inertial 
frequency was shifted about 3-6% higher than the inertial 
frequency and it dominated the CW spectra calculated fron 
the mixed layer velocity measurements at CMA3 (Fig. 14). 
The rotary spectrum estimates for “he inertial period 
moticns are given in Table IV. The CW, inertial HKE vari- 
ability estimates are generally 1 to 2 orders of magnitude 
larger than the CCW rotating motion. In comparison, the 
subinertial (pericds of 3 IP) and superinertial (periods of 
12 h) CW spectral estimates ars more than an order of magni- 
tude less than the inertial moticn (Tables VII and VIII). 

At all levels at the deep region, the positive rotary 
coefficients approached unity, which indicates that the 
waves were polarized in the CW direction (Table IV). The 
direction of tte semi-major axes of the inertial ellipses at 


CMA2 rotated CW with depth except near the botton. The 
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TABLE IV 
Normalized Rctary Spectrum Analysis at Inertial Frequency 


Normalized 
Rotary 
Spectra Ellipse 
Lepth CW CCW Der. Rotary 
(n) (cm/s) 2/cph Stab. (Deg.) Coeff. 
a) CHAI 
49 6.7x10*3 1.3x10+3 0.81 18 *0.68 
64 2.5x10+3 1.4x10+3 0.20 65 +0.62 
b) CHAZ 
19 1.8x10*5 3.5x10*3 0.90 60 *0.96 
179 4.6x10++ 3. 7x10+2 0.95 80 +0.95 
224 2.8x 10++ 1.3x10+2 0.30 60 +0.96 
c) CNA3 
21 4.0x10#5 2.0x10*3 0.87 50 +0.98 
251 1.4x10+5 6. 1x10+2 0.63 90 +0.98 
437 1.2x10*5 6. 3x10+2 0.50 150 +0.98 
457 6.2x10** 52.0x10:“ 0.77 160 +1.00 


direction of the near-bottom ellipse rotated CCW with depth 
between the +termcciine (179 m) and near th= botton (324 n). 
The stability cf the ellipses indicated that the inertial 
moticn near tke bottom at CMA2 was isotropic. In cCOHUE£€”, 
inertial cscillations were strongly anisotropic at all of 
the cther depths in the deep region. The inertial vave 
response at CMA1 was markedly less intense than at CMA2 or 
SUAS, but it was polarized ina CW sense. The reasons for 
the semi-diurnal tide at CMA1 being much more energetic than 


the inertial nction is explained below. 
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The vertical and horizontal coherences are computed to 
estimate the spatial scales. The spatial scale L is con- 


puted from the cross-spectral phase relationship using 


l= == fe ә (3) 


where Ad is either the vertical or horizontal distance 
cen the current meters and Ad? is the phase difference 
between siailar velocity ccmponents (Pollard, 1980). This 
analysis is restricted to the CMA2 and 3 due to the absence 
Of an energetic inertial response at CMA1 (Table V). Esti- 
mates of vertical scales are more accurate than the horizon- 
tal scale estimates because the vertical spacing of the 
current meters below the sixed layer is not the same at the 
two arrays. These scales calculated from the entire storm 
time series will be compared belcw with the scales calcu- 
lated from the instantanecus phases. 

The vertical cross-spectral analyses (Table V) indicate 
that the inertial response in the mixed layer led the nsar- 
bottom response. The normalized cross-spectrum variance 
between the vertical levels was about a factor of two 
greater at CMA3 than at CM22. The vertical scales computed 


from (3) are about 360 and 540 m at CMA2,3 Tespectively. 
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That is, the inertial motion had scales of the order of the 
water depths at both arrays and was coherent at the 95 2 


confidence level. 


TABLE V 


Vertical and Horizontal Coherencies at the Inertial 
Frequency 


Normalized 
Cross-Spectra 77: Phase 
Array Variatles Varianca/cıh Coherence (deg.) 
vertical 
2 u(s) ,u(b) A 0.62 -50 
2 О ли 2, 710 ** 0.60+ -70 
3 u (s) ,u(b OSS XU 0.70 O 
3 v (s) ,v(b) 6.Cx10+4 0.80 -50 
Horizental 
5,2 u (s) 1. 9x10“5 0.98 -45 
3,2 ç (S 1.6x10*5 1.00 -55 
5,2 Y ib) 2:5x1077 0.98 -20 
72 v (b 2.1x10+* 0.95 -30 
memes Significant at 95 % confidence 
* spectra assumes an average .. frequency of 0.0412 cph 
| = surface layer (A = 19 461600 21 n] 
İb = bottcm layer (CMAZ = 324° m ; CMA3 = 457 nm) 


The instantaneous phase differences of the velocity con- 
ponents between the two records can also be used to estimate 
the vertical scales (Fig. 15) Initially, the estimated 
vertical scales at CMA3 were much greater (typically 800 to 
900 m) than the water depth, but after UD 257 the estimates 


were consistently of the order of the water depth. The 
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large vertical scales observed initially agree with thosa 
estimated from the Eloise simulations by Price (1983) and 
the Allen cbservations by Erooks (1983). At CMA2, the ver- 
tical scales had similar trends and converge to scales of 
the crder of the water depth. 

The horizcntal cress-sfectrum estimates (Table V) indi- 
cate both velccity components were coherent in the mixed 
layer. The corresponding scales, as estimated from (3), are 
20 and 35 km inthe east-west and north-south directions 
using an array separation cf 34 km. Initially, the horizon- 
tal scales were set by the atmospheric forcing, but these 
estimates may be somewhat biased because the entire stcrm 
time series was used in the computation of che spectra. The 
cross-spectrum variance between the two near-bottom records 
was about cns order of magnitude less than in the mixed 
layer. However, these values are significant at the 95 % 
confidence level as indicated by the coherence squared. The 
respcnse at CMA2 leads the response at CMA3 by about 2 h. 
The horizcntal scales near the bctton calculated from (3) 
are 20 and 35 km in the east-west and north-south direc- 
tions, respectively. The differences in the scales are due 


to the larger scales set up by the hurricane in the 
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Eong--rack (ncrth-south) than ne Cross-tzack  (east- 
west) directions. These estimates of scales may be somewha- 
misleading because of the instrument depth difference (135 
m), although tkey agree well with mixed layer estimates. 

The instantaneous phase differences between the mixed 
layer currents indicate that the horizontal length scales 
were set by the storm initially (Fig. 16). The scales 
decreased rapidly and approached constant values of 60 and 
100 km in the east-west and north-south directions, respec- 
tively. These scales in the mixed layer are roughly 1-2 
times the karoclinic Rossby radius of deformation for the 
first mode (56 km). The estimates of the horizontal scales 
of inertic-gravity waves in the tharmocline are computed 
from the tottci and mid-depth instantaneous phase differ- 
ences in the U-velocity ccmponent at CMA2 and CMA3, respec- 
2 “ely (Fig. 17). The east-west and north-south scales in 
the thermocline immediately after hurricane passage are 
roughly the same as in the mixed layer. After 4 IP, the 
scales revert *c the pre-storm scales of 25 and 50 km in the 
east-west and north-south directions, respectively. Scme 
cauticn is advised in accepting these estimates due to the 


large depth difference of 73 m. Taa later section,- it is 
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suggested that the energy was contained primarily in the 
barotropic and first two taroclinic modes. Hence, a depth 
difference of 73 m may be acceotable wi-hin the dəeper ther- 
mocline because of the relatively large displacement of the 
isotherms, tre large vertical scales, and the amount of 
energy contained in the barotropic mode. 

The estimates of vertical energy propagation speed are 
computed only cn the basis of the lag times between the sur- 
face and suksurface response (Table VI). In the first case, 
the response is felt at the bottom of the water column 
Within 3 to 4.5 h after increases inthe mixed layer are 
detected. Tte initial respense corresponds to a vertical 
scale of roughly 900 n (Fig. 15). This lag 0£73 to 15574 
results in a vertical energy prcpagation speed of 2.6 kn/d 
which is rcughly 1 to 2 orders of magnitude greater than the 
values calculated by Brooks (1983) and Price (1983). The 10 
IP lag in the second maximum observed at 251 m represents a 
vertical enercy propagaticn speed of 26 m/d when the verti- 
cal scale is about 500 m. This vertical propagation of 
energy agrees quite well with the estimates derived from the 
Allen observations. However, there is also a second maximum 


observed near the bottom about 6 IP following the storm. 
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Consequently, the vertical (Cga) and kOrizontal (Сан ) 


group velocities are calculated from the following expres- 


sions (Brocks, 1983) 


2 
" G N 2 
— ә 2 tan 8 9 (4) 
2 
= O N 2 E 
Cg = k oz tam — 7 (8) 
where o is the cbserved frequency, m and k are vavenumbers 


m” she vertical and horizcntal cross-track directions, and 
tanƏ = k/m. The wavenumbers m and k are equal to 27 /L4,and 
gir /Ly respectively, where L, and L, represent the vertical 
and horizontal scales as estimated previously. The assigned 


values to the variables in (4) and (5) are 


Q 
II 
J 

e 
Un 
= 
” 


10-5 rad/s, 
N2 = 3.30 x 1C-5 rad?/s?, 


п = 6.98 x 10-3 rad/n, 


Il 
-— 
° 
=à 
O 


m x 10-2 rad/m, 
k = 1.05 x 10-* rad/m. 
The initial vertical group velocity is 1.25 ka/d as com- 


puted from tte dispersion relation (4) wi- han initial 


scale of 900 m (Table VI). This estimate is about half of 
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TABLE VI 


Ccmparison cf Vertical and Horizontal Group Velocity. 





+ 


Vertical Brooks * 
Scales (1983) 
5600 m 


— ng 
—11 
Wr: 
С0 О 
Cə 
us 


900 m 
Vertical 
Grcup, 
Velocity 


Lag Tine 2.6 0 
Dispersion 1.3 
(kn/d) 

Berizontal 

Grcup. 

Velccxty 


Dispersion 80. 30; 23% 86. 
(km/d) 


| 
* From Brooks (1983) | 


3 
7 “06 26: 


— — — “ . t J. w Qasa s ¿sam “m sme Q ei —= Xə 


the cne ccmruted based on the initial lag tine and is still 
an order of magnitude larger than the estimates computed 
from the Eloise simulations. After a few LP, "the vertical 
S CUE velccity is about 0.27 km/d, which is of the order of 
the Eloise estimates cf abcut 0.16 km/d (Price, 1983). Fur- 
thermore, this estimate is about one order of magnitude 
larger than the value of 0.02 km/d that Brooks (1983) esti- 
mated from tke Allen observations. Clearly, the initial 
Storm forcing causes energy to propagate vertically at a 
Significantly Eigher rate than after a few IP following the 


- ora. 
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The  hcrizon*tal propagation of energy is a much more 
Sapa than the vertical prcpagation (Table VI). The grcup 
velocity estimate for the Frederic observations is initially 
80 km/d; and it decreases to abcut 30 km/d a few IP follcw- 
ing the stcrn. This difference is due to the tan? depen- 
dence on m as previously defined. Tho “initial Neri zontal 
group velocity estimate agrees with the Price (1983) esti- 
mate whereas the lower value agrees with the Brooks (1983) 
estimate. Trerherizcntal phase speeds of the inertio-grav- 
ity waves in the wake of Frederic are about 70 cm/s with a 
vertical phase speed cf rcughly 1 cn/s. 

In summary, inertial waves maintained a constant ampli- 
tude over 1 to 2 IF fcllowing the passage of Frederic. 
These oscillaticns in the mixed layer then relaxed over an 
e- folding scale of about 4 IP. In the subsurface layers, 
the forced inertio-gravity waves decay2d over similar time 
scales after tke appearance of a secondary maximum. The lag 
between the initial impulse of energy and the second subsur- 
face maximum varied over depth. For example, this period 
was about 10 IE at 251 m, but it was only 6 IP at 457 n. 
These lag values define a modulaticn envelope which evolves 


in time. If vertical prepagation of energy from the mixed 
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layer alone acccunts for these secondary maxima, why do the 
maxima occur later in the deep thermocline than near the 
bottcm ? Cther mechanisms besides dispersion of energy from 
the mixed layer must be responsible. For example, nonlinear 
resonant interactions among the  inertio-gravity waves could 
be responsible (McComas and Bretherton , 1977). Analyses of 
these interactions are beycnd the scope of the present work. 

The initial vertical scales of the inertio-gravity waves 
in the wake of Frederic are approximately equal to the 
Eloise and Allen estimates. Hcwever, after a few IP the 
vertical scale decreased to about 500 m or approximately the 
water depth. The horizontal scales estimated from the Fred- 
əric cbservaticns differ qreatly from the Allen data and 
Eloise simulations. The hcrizontal scales during Allen were 
four times the mixed layer horizontal scales observed inmed- 
iately after the passage of Frederic. This difference 
increased tc a factcr of about eight approximately 4 ID 
after Frederic. These discrepancies may be due to the rug- 


ged Ecttom topography in the DeSoto Canyon region. 


B. SUBINERTIAL RESPONSE 
The forcing ty hurricare Frederic also altered the "mean 


flow" and the lower frequency motions (periods greater than 
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3 IP). The variability and respcnse of the mean flow to the 
forcing is depicted in progressive vector diagrans (PVD). 
The ccean current data were low-pass filtered with a half 
power point at 26 h to isolate the forced response on the 
lower frequency variability. 

The mean flow in the mixed layer at CMA2 was directed 
towards «he east pricr to the stcrm (Fig. 18a). During the 
passage of hurricane Frederic near JD 256, there was only a 
small deflecticn in the mean flow. The subsequent horizon- 
tal displacements indicate that the mean flow was about 26 
km/d, or almost twice the mean current during the pre-stcrn 
pericd. The mean flow at 179 m (not shown) at CMA2 also 
increased towards the east. 

The bcttcm PVD at CMA2 indicates that the mean flow at 
that level was mcre or less constant (Fig. 18b). However, 
the direction of this flow was opposite to that observed in 
the mixed layer. Consequently, the bottom PVD was nearly a 
nirrcr image of the surface PVD, although the bottom flow 
was slower than that cf tke mixed layer. There was also a 
superposed lower frequency motion of about 3 IP near the 
bottom, particularly durirg the period between the initial 


and secondary maxima of HKE (first 6 IP in Fig. 9). The 
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lower frequency oscillaticns may be due to coastally trapped 
waves, which were enhanced by the storm, since waves of the 
same period were observed during the quisscent period 
(Appendix E). These waves can be altered by both botton 
topography and continuous density stratification (Wang and 
Mooers, 1976). A study cf the mechanisms by which the mean 
flow is altered by these waves and the rugged terrain is 
also keyond tke scope of this thesis. 

At CMA3, the surface mean flow was southward pricr to 
the storm (Fig. 19a). At the time of storm passage, near JD 
256, the surface mean flow was deflected outward from the 
storm center. Between JD 257 to 269, the mean flow acceler- 
ated to about 13 km/d and then decelerated to about 6 km/d 
towards the scuth. Prior to the storm, the mean flow near 
the Ecttom at CMA3 was ncrthward following the local bottom 
topography (Fig. 19b). There was an eastward (up the canyon 
slope) deflection late on JD 256. After completing a loop, 
the mean flow veered towards the northwest (down the canyon 
slope). Between JD 254 to 263, the mear flow near the bct- 
tom was considerably faster than either earlier or later 
times. Again, the mean flow near the botton vas reversed 


relative tc that in the mixed layer. 
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Scale: k 
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EC Scale: A 
70 km 


Figure 18. Lcw-pass Filtered PVC at CMA2. The * represents 
a time interval of 48 h 2:5 on JD 248 for 
R mixed layer (19 n) and D) -tom Layer (324 
m). 

Thermocline temperatures in the wake of Frederic varied 
spatially as well temporally. There were sizeable subiner- 
tial oscillaticns as well as near-inertial oscillations in 
the temperature signal (Fig. 20). These temperature varia- 


tions were typically 3.09C at both array sites and were 


superposed on mean temperatures that range between 13.0 to 
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Figure 19. Lcw-rass Filtered P C 
tarsıing on JD 285 for 


VD a 
a time interval of 48 h 
A nixed layer (21m) a 


t 
s 
nd b) bottom Layer (457 


14.59C in the thermocline. At CMA3, the temperature at 251 
m started to increase about 0600 GMT on JD 255 and continued 
+o increase over the next 15 h due to the strong northward 
current. The temperatures then decreased by 29C due to a 
strong southward current. Ineremal fluctuations of abcut 


39C persisted cver the subsequent 8 IP. After this time, 
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the variations undervent a slight increase in frequency 
which corresgcnded to about the time of the secondary maxi- 
num in the horizcntal velccity components. For the remain- 
der c£ the record, the temperature variations were 
superposed on a staticnary mean. 

Ecttom temperature variability at CMA2 is a clear exam- 
ple cf advective processes that are associated with lower 
frequency moticns (Fig. 21). On JD 254 (11 September), a 
temperature decrease of 1.39C over a 16 h period marked the 
beginning cf the fcrcing period as cool water was advected 
from the DeSoto Canyon. The temperature then increased by 
3.79 C over tke next 25 h as warm water was advected by the 


southward velocity component which was downslope in this 


part Of the canyon. As the winds relaxed, cooler water was 
again advected into the area. Around JD 260, the tenpera- 
ture dropped ty 49°C cver an 11 hour period. Tus drastic 


change in temperature was due to the advection of cold water 
by these subinertial waves of 3 IP period. From the orien- 
taticn of the ellipses (Table VII), the direction of the 
wave propagation was northward (upslope, away from the 
DeSotc Canyon). These waves set up cross-shelf oscillaticns 


(relative to CMAJ 5289 perno Gth-scuth direction. After 
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these large variations,  near-inertial oscillations in the 


temperature were superposed on a gradual warming trend. 


TABLE VII 
Normalized Rotary Spectrum Analysis at Subinertial Frequency 


Normalized 
Rotary 
Spectra Ellipse . 
Depth CW CCW Dake. Rotary 
(m) (cm/s)2/cph Stab. (Deg.) Coeff. 
324 4.5x10+3 3.041073: 0.65 19 +0.20 
457 4.7x10*3 1.0x10+2 0.50 10 +0.98 


The potential energy spectra ware computed using the 
+amperature time series and the vertical temperature gradi- 
ents from the AXPT data ccllected by Black (1983) (Appendix 
A). The mcs+ prcminent peak in the potential energy spectra 
at CMA3 is near the inertial frequency (Fig. 22). The most 
energetic peaks at the semi-diurnal and subinertial frequen- 
cies are near the botton. Potential energy variability in 
the bottom layers at CMA2 (not shown) is significant at the 
95 % confidence levels only for periods of about 3 IP. 

From rctary spectrum analyses at periods of 3 IP near 
75 Ectton: 177 ::6€ CF and CCW rotating components at CMA2 
were nearly equal, with just a slight tendency towards CW 


polarization; and (2) the motion at CMA3 was polarized ina 
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Figure 22. 
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CW sense (Table VII). Further more, the phase propagaticn 
for these lcnger pericd waves was northward, and anisotropic 
as indicated Ly the ellipse orientation and stability. 

The cross-spectrum between the temperature ard north- 
south velccity records at the bottom was also dominated by 
the peak at about 3 IP (Fig. 23). The corresponding phase 
and ccherence estimates are 100 deg. and 0.98 respectively, 
with the velccity leading the temperature signal. Hence, 
these lower frequency oscillations in the cross-shelf direc- 
tion at CMA2 were properly oriented for advecting the temp- 
erature pattern. These cross-shelf oscillations verə 
respcnsible fcr the variability in the bottom temperature 4 


IP after the passage cf Frederic. 


e SUPERINERTIAI RESEONSE 

As Frederic passed within 130 km to the west of the 
CMA1, the dominant  respcnse was a significan: increase in 
the westward current (Fig. 24). This increase was due to 
convergence of flow on the right side of the storm. By ccn- 
trast, there was a divergence of flow on the left side of 
the hurricane, where water was transported off-shelf by the 
Winds. For example, a negative tide was recorded at Biloxi, 


Mississippi which was indicative of the left side of the 
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stormn. The north-south current only increased by about 20 
cn/s, whereas the east-west current increased to 80 cmys. 
After the hurricane passage, near-inertial waves persisted 
over the subseguent 3 IP. The post-storm semi-diurnal tidal 
currents vere mcre energetic than during the quiescent 
pericd and were superposed on a non-stationary trend. At 
the intermediate depth at CMA1, the current speed increased 
B 0 CM/S during the period of strong forcing, as the 
semi-diurnal tidal currents were enhanced. 

The presence of internal tides at the semi-diurnal fre- 
quency is indicated by the HKE spectra during the storm 
pericd a+ CMA1 (Table VIII). The spectral estimates in the 
semi-diurnal frequency band are almost equal *o the inertial 
p2ricd estimates. Baines (1973) shovs that one necessary 
Eordition for tke 5255 of internal tides is a marked 
increase in the semi-diurnal tidal currents. The observed 
increase in tte semi-diurnal tidal currents should be equal 
to or greater than the inertial effects. The storm-spectral 
estimates at the semi-diurnal frequency exceed those from 
the quiescent period by mcre than order of magnitude (Appen- 
dest B). A barotropic tide propagating over the shelf break 


generates internal tides cf semi-diurnal tidal period 
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Temperature from CMA2 at 324 m. K, represents 
the diurnal tide frequency, f the inertial 
frequency and M, the semi-diurnal tide 
frequency. 
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Figure 24. Mid-depth Velocity Components from CMAI. The 
ordinate depicts magnitude in cm/s of “he east- 
west currents (upper) and north-south currents 


(lower). 
(Prinsenberg and Rattray, 1975 ; Barbes et al., 1975 ;  Tor- 
grimson and Hickey, 1979). Hence, an increased barotropic 
tide associated with the storm surge propagating over rugged 
bottcm topography should enhance the internal tides as nani- 


fested in the semi-diurnal tidal currents. 
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ТАВИЕРУТ11 


Nermalized HKE 5. at Inertial and Semi-diurnal \ 


Frequency 
Normalized 
Meter NO. Of Data Fr equency su — 
Lepth Points Resclution 
(m) (cph) (EnZs) 27eph 
CMAN 
49 3200 0.0018 2.0x10*3 1509 10:83 
64 2410 0.0025 3.0x1017 6.0x10+2 
92% 3000 0.0020 1.4x10+3 32201074 


* instrumentation problems with time clock 
ae the semi-diurnal tide frequency band 
f is the inertial/diurnal tide requency band 
The dcminance of internal tides on the shelf at CMA1 can 
be explained by considering the bcttom slope in the region. 
If the bottcm slcpe is greater than the internal wave char- 
acteristic, then the rays will not be admitted to the shelf, 
but instead will be reflected into the oceans interior (Bar- 
beleWe- al, 1975 ; LeBlond and Mysak,1978 ; Torgrimson and 
Hickey, 1979). The critical slope for the internal wave 


characteristic is given by: 


ca = ) Co) 


where f is the local Coriclis parameter, which is equal to 


7.18 x 1075 rad/s near the DeSotc Canyon. The remainder of 


77 






the variables have been defined above. The slopen of the 
inertio-gravity wave characteristic is about 4 x 10-3, 
whereas in the direction of the canyon axis (45 deq.) the 
bottom slope is 6 x 10-3 . Therefore, it is clear that the 
inertio-gravity waves will be reflected seaward as the slope 
of the characteristic is less than the bottom slope in the 
DeSctc Canyon region. Inertial oscillations will be locally 
generated on tke shelf but they will only persis* for a few 
IP fcllowing tbe storm. Ihe internal tides, however, will 
be admitted to the shelf region because “hə slope of the 
semi-diurnal tide characteristic is 2 x 10-2, which is 
almost an crder of magnitude larger than the bottom slops. 
The semi-diurnal *ides were also enhanced at CMA2 and 3, 
but they were much less energetic than the inertial waves 
(Cable IX). The CW polarized, semi-diurnal tidal variabil- 
ity was more than an crder of magnitude greater than the CCW 
variability estimates, except in the bottom layers a+ CMA3. 
The directicns of the semi-major axes of these waves changed 
from 120 tc 20 deg. with depth at CMA2, whereas the opposite 
occurred at CMA3, where the direction changed from 30 to 50 
deg. Hence, the waves were anisotropic and exhibited a 
E ng north-south cemponent of flow at CMA2 and to the 


northeast at CMA3. 
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TABLE IX 
Normalized Rctary Spectra at the Semi-diurnal Frequency. 


Normalized 
Rotary | 
Spectra Ellipse 
Depth W CCW Dir, Rotary 
(m) (cm/s) 2/cph Stab. (Deg.) Coeff. 
a) CMA2 
19 1.3x10+3 4.3x10+t2 0.80 120 6010 
179 5.0x10+2 O 60 +0.72 
324 4.3x10+2 Wzxı 082 0.87 20 +0.41 
Ы) CMA3 
21 1.0210*> 5.5x10+2 0.58 30 +0,55 
251 6.7x10+2 1.0x10+2 0.92 90 10:70 
437 3.2x10+2 1.3x10+2 0.40 40 +0.40 
457 2.0x10*? mer 102 0.15 50 +0.04 


In summary, the superinertial responsa in the wake of 
Frederic was c¢cminated by the enhanced sami-diurnal tidal 
current at all three arrays. The scemi=diúrnàl tidal cuz- 
rents were markedly more energetic than the inertial cur- 
rents on the continental shelf at CMA1, and agrees with 
internal tide thecry (Baines, 1973). The storm surge, which 
is associated with the passage f the storm, increases the 
BXrotropic tide. As this increased barotropic tide propa- 
gates over rucged bottom terrain in a stratified ocean, the 
increase in tte semi-diurnal tidal currents are manifesta- 
tions of internal tides (Barbee tal, 1975 ; Prinsenberg 


w Rattray, 1975 ; Torgrinson and Hickey , 1979). Internal 
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tides can propagate freely up the continental slope and onto 
the shelf, whereas inertio-gravity waves generated off the 
slope are reflected tcwards the interior of the ocean. The 
rotational characteristics of the semi-diurnal tidss, as 
well as the direction of phase propagation of these waves, 


are consistent with these theories. 
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A. THEORY 
The Sturm-Licuville pıcblem is solved using the vertical 
profile of Brunt-Vaisala frequencies and linear wave theory. 
The vertical eigenvalues cf the horizontal velocity are then 
used to oktain a least squares fit to the Fourier coeffi- 
cient time series computed from the demodulation of the cur- 
rent meter records. The amount of variability associated 
with each mode is detərmined to assist in understanding the 
modal structure of the inertio-gravity wava response gener- 
ated ty the passage of a hurricane. 
The problem is formulated by following the work of 
Fjeldstad (1958) and making the following assumptions: 
e f plane; 
w continuous stratification; 
e hydrostatic balance with basic state at rest; 
ә jincomrressikle; 
e inviscid; and 


e flat kottcn. 
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The governing vertical structure equation for linear, iner- 


*io-qravity waves is giver by 


d^ 


dz 


-— 





* u^o(z)W - 0 (7) 


DI 


where W(z) gives the vertical structure of the inertio-grav- 
. Í 3 x 

lty wave, Riz K7o?-£ e $(2)5N*- c2, and Ke = ke + 12, and K 
and l are the horizontal wave number vectors in the x and y 
directions, respectively. Furthermore, equation (7) is a 
seccnd order,  hcmogeneous, ordinary differential equation. 
The specification of the problem is completed by imposing 


the followinq boundary conditions: 


dW » = 
dir guW = 0 2=0 9 (3) 
W = 0 z--d ə (9) 
The surface Ecundary conditicn (8) is a dynamic boundary 


condition which states that at the sea surface pressure is 
continuous across the interface. Boundary condition (9) 
specifies ne normal flow through the botton. The solution 
to eguations (7-9), which define a well known Sturm-Liou- 
ville problem, yields the structure of the vertical velocity 
for ccnstant N, and, is given by 


4 nuz 
W (z) 5-52 3 ) 2 010) 
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where the mode number n=0,1,2..., The horizontal velocity 
eigenfunctions U (z) are given by 


U (2) x cos( E) ? E) 


which is the vertical derivative of W, (z). 

The Sturm-Licuville problem is solved numerically by the 
predictor-corrector fourth-order method modified by Hamming 
(Gerald, 1980). This methed computes a new vector from four 
preceding values. A fourth-order Runge-Kutta method is used 
for the adjustment of the initial vertical increments and 
the ccnputaticns of starting values. 

The resulting amplitudes of the horizontal velocity 
eigenfunctions fcr modes 0 through 3 are shown in Fig. 25. 
The number of zero crossings equals the mode number. Mest 
of the vertical structure lies in the upper thermocline 


regicn near the tase cf tte mixed layer. 


B. LATA ANALYSIS 

The horizental eigenfunctions are then fitted in a least 
Squares sense tc the time-varying amplitudes of the horizon- 
tal velocities. The mathematical details are given in 
Appendix A with the final matrix equations (16-18). The 


time evoluticn cf the north-south velocity coefficients at 
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CMA3 indicates that the storm excited moies 0,1,and 3 witha 
damped mode 2 (Fig. 26). The amplituies of modes 0 and 1 
reached a peak on JD 256 and then decreased. By JD 262, the 
barotropic mode horizontal speed decreased to nearly a con- 
Stant value of about 10 сп/5. After JD 262, modes 0 and 1 
oscillated over the remainder of the record with a period 
of about 5 IP. This behavior of the mođəs in Fig. 26 indi- 
cates that a modulaticn occurred between in=rtial wave modes 
0 and 1. More importantly, these modes define a modulation 
envelope which is consistent with the amplitude variations 
of the currents below the mixed layer. 

At CMA2, the time-varying coefficients of the modes are 
markedly different in their behavior (Fig. 27). The spatial 
variability in the modal structure between CMA2 and CMA3 is 
associated with the differences in the mean flow, which is 
topographically controlled. Initially, all of the modes 
were excited, with a fairly energetic node 3. The baro- 
tropic coefficient decreased to about 2 cm/s while the mode 
1 amplitude approached zero near JD 252. Arteon JD 262; 
modes 0 and. 1 were in phase and defined the linits of the 


coefficients. 
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Bada current time series are re-computed using the 
coefficients from the least squares fit and the eigenfunc- 
tions for each cf the current meter depths at CMA3. Eos 
example, the reconstructed time series are compared in Fig. 
28 with the actual demodulated current time series at 251 m. 
Recall that a seccndary maximum occurred at this site about 
10 IE follcwing the passage of Frederic (Fig. 10). The 
reconstructed amplitude cf the barotropic mode is greater 
than the actual v component of the data between JD 255 and 
moe (Fig. 28). However, Summing modes 0,1 and 2 accounts 
for most cf tke cbserved variability, with the barotropic 
mode the most energetic. Over the remainder of the record, 
80 0 tc 270, most of the observed variabili-y can be 
described cnly Ey modes 0 and 2. Diewscontr.bucerson of the 
barotropic mode was significant and consistent with earlier 
estimates that the vertical scales were of the order of the 
water depth. The strength of the barotropic mode continued 
for the entire time seriss of the modal coefficients at 
CMA3. 

Estimates of the modal contributions to the observed 
near-inertial variability are given in Table X for both the 
storm and poststorm periods at CMA2,3. These estimates are 


based on the expression 
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AR (12) 


where Y; and Vj are the observed and modal scalar components 
25 ne horizcntal velccity. The index m depicts the summa- 
5002 of the mcdes, for иә, 0, 0+1, anà 0+1+2, and index 
j repesents a summation cver the number of observations in 
the period. The storm period is redefined to start at the 
time of hurricans passage and continue for 7 IP. The post- 
storm period starts where the storm period ends and contin- 
27 7 oC rouqhly 8 IP until the end of record. 

The mixed layer variability at CMA2 is dominated by the 
barotropic mode during the storm period. The contribution 
of the barctrcpic mode tc the cbserved horizontal current 
variability decreases with depth, but by including the 
baroclinic mcdes most of the near-inertial variance can be 
described ty modes 0,1 and 2. The only exception ís at 179 
n. The first three modes contribute only 54 and 37% to the 
observed hcrizcntal current variance in the east-west and 
north-south directions, respectively. Furthermore, the 


u-velccity  ccmpcnent of the barctropic mode exceeds the 
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TABLE X 
Variance cf the Ncrmal Modes at CMA2 and CMA3. 


o 021 0+1+2 
Depth Velccit +1 +1+ 
() Çonponeht Bericd (2) (2) ( 
CMAZ2 
19 u s 68 78 100 
19 u ps 87 91 100 
19 V 69 78 100 
19 V ps 22 93 100 
179 u S ə. 50 54 
179 u ps 98 98 95 
179 V S 35 31 37 
179 V PS 97 98 95 
324 u s * 50 67 
324 u ps 97 98 99 
324 V S 4 5 72 85 
324 V ps 97 99 99 
CMA3 
21 u s 69 97 100 
21 u ps 73 86 100 
21 V S 68 98 100 
21 V ps = 50 100 
251 u s 55 81 98 
251 u ps 95 96 29 
25] V s 58 8ü 98 
251 V PS 94 95 99 
437 u S 73 68 97 
4:7 u ps 95 94 91 
437 V s 64 78 94 
4 37 V ps 94 94 92 
457 u s 22 82 85 
457. u ps 79 74 99 
457 V S 21 89 84 
457 V ps 81 66 99 
s: storm pericd ( number of otservations = 1008) 
ps: post-stcrm pericd (number of observations = 1152) 
*: an overestimation of the variance (explained below) 


observations by a considerable amount and causes equation 
(12) to be less than zerc (see Fig. 227. This behavior 
could be attributed to the botton boundary layer vhere iafi- 
nite energies occur when the slope of the internal wave 
characteristic approaches the bottom slope (Prinsenberg and 


Rattray, 1975). Most of the observed variance during the 
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post-storm period is attributed to the barotropic mode. The 
modes are much more well-tehaved during this period, which 
is presumatly a manifestation of free waves. 

The modes of the horizontal velocity are well-behaved 
Being the stcrm period at CMA3, except near the bottom 
where only 20 to 22% of the observed variance can be attrib- 
uted to the tarctropic mcde. This is probably due to bottom 
boundary effects. At the other depths, the barotropic mode 
plains 55 <c 73% of the variance, with most of the vari- 
ance accounted for by the summation of the barotzsopic and 
Gist two tarcclinic modes, During the post-storm period, 
the cbservaticns of the mixed layer velocity are less than 


Ehat cf tke tarctropic mcde in the north-south direction, 


1 


vhich causes the value to te less than zero. The reason £o: 


t 


iD 


this large discrepancy is due to the e-folding scale of th 
horizcntal velccity in the mixed layer, the persistence of 
20 ccefficiənt cf the barotropic mode (see Fig. 26), and 
the radiation of inertial waves away from the storm track. 
Otherwise, 72 tc 95% cf tte observed horizontal current var- 
iance can be attributed to the barotropic mode, and adding 
the first twc taroclinic mcdes accounts for more than 90% of 


the variance. 
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Hopkins (1982) modeled only the baroclinic response <o 
hurricane passage and he found that the ELE 
waves decayed toc rapidly in the thermocline and bottom lay- 
ers to account for the vertical structure of the observed 
variability. Thus, a model with only the baroclinic modes 
could not explain the secondary maximum observed in these 
layers, 

Some caution h2s +9 applied in the interpretation of the 
normal modes for several reasons. First, tne formulation of 
the Sturm-Liouville oroblem assumes that the vertical struc- 
ture consists of a standing wave for the N2 = constant case. 
rrom kasic physics, a standing wave caa be decomposed into 
two waves of equal amplitudes propagating in opposita direc- 
tions, which assumes no phase propagation. Over the first 7 
to 10 IP, there is a downward propagation of energy, as well 
as upward propagation of phase, which violates the standing 
wave argument. Secondly, a flat bottom ocean is assumed for 
the application of normal mode theory. The validity of this 
assumption crucially depends on the ratio between the slope 
of the internal wave characteristic to the bottom slope. 
Por example, as this lratic approaches unity, the energy den- 


Sities in the bottom boundary layer approaches infinity. At 


92 





245 250 299 260 265 270 


20 


MODE O 


V—Velocity 
(cm/s) 
= 


V-Velocity 





0 
20 
m 
ар ти M... 
= MODE 0+1 
0 
20 
= ! 
> 2.70 
Sm 
=~ 10 | 
p š MODE 0+1+2 
> — 
9 e 
245 250 255 260 265 270 


Time (Julian Days) 


Series at CMA3 for the observed 
d) and the reconstructed time 
O 


time series ) 
r a) node 0 , b) mode 0:1 , c) 


series 75 
mode 04 D 


Figure 28. V Component Tine 
(scli 
ed) f 


93 





sat point, dissipation becomes important and must ba 
included in the medel (Prinsenberg and Rattray, 1975). Con- 
sequently, the assumption of inviscid dynamics is no longer 
valid. Despite all these approximations, the normal mode 
analysis fits the observations quite well. Even the reso- 
nance or modulation envelcşe which occurred at CMA3 about 10 
IP fcllowing tke storm passage, is represented quite well by 


modes 0,1 and 2. 
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VI. CONCLUSICNS 


mr. 


The ocean response near the DeSoto Canyon to hurricane 
Frederic was dcminated by the excitation of inertial waves 
in the mixed layer and inertio-gravity waves in the thermo- 
line. These waves are not admitted onto the shelf because 
*he slope cf tte characteristics cf internal wave motion is 
less than the bottom slope (Barbee et. al, 1975 ; LeBlond 
and Mysak, 1578 ; Torgrinson and Hickey, 1979). The major 
features cf the forced inertio-gravity waves observed at 
25052 and 3 are: 


e initially, the vertical propagation of “57 was of 
the order cf 1. km/d as the forcing was felt throughout 
the water column within 3 to 4.5 following the pas- 
sage cf Frederic; 


e the initial horizontal propagation of energy was abcut 
80 kmyd with a phase speed of 70 cm/s; 


ə the waves were anisctropic at all depths and had an 
upward propagaticn of phase; 


° maximum HKE in the mixed layer occurred at the fringe 
cf the maximum wind regime with a second maximum in the 
thermocline 6 to 10 IF after the storm; 


e che EKE decayed on an e-folding scale of about 4 IP 
after the maximum in the mixed layer and the seccnd 
maximum in the thermccline; 


Scie Vertical propagation after about 6 IP is of the 
order of 102 m/d with a corresponding horizontal propa- 
gation of energy of 30 km/d; 


e the mixed layer horizontal scales were 60 and 100 km in 
the cress and  along-track directions, 7::5-- 
while the scales in the thermocline were 25 and 50 kn 
in. he CrCSS and along-track direceions, 
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e -the vertical scales were cf the order of the. water 
depth which suggests that the barotropic mode was 
inportant; 


e most cf the observed variabilty was in the barotropic 
mode, with some, ccntributions from modes 1 and 2; 
together they defined a modulation envelope in the deep 
thermecline at CMA3; 


e the modal ccefficients varied in space as well as time 
due tc tte rugged bottcm topography and mean currents. 


The mean flcw was topcgrapically controlled and reversed 
w ection frcm surfaca tc botton. The forced response in 
the mixed layer caused the mean flow to increase at both 
CMA2 and CMA3. The mean flow below the mixed layer was also 
changed by «the storm. Furthermore, this mean flow was 
influenced ty subinertial variability (periods of 3 IP) and 
was linked to the advecticn of cooler water cross-shelf rel- 
Eye “o CVA2. 

The semi-diurnal tidal currents increased in magnitude 
and varied spatially in response to the passage of hurricane 
Frederic. The near dcminance of the semi-diurnal tidal cur- 
rents over the inertial mctions at CMA1 was consistent with 
internal tide thecries. The semi-diurnal tides were admit- 
ted to the shelf region because the slope of the internal 


wave characteristic exceeded the critical bottom sicpe 
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fSa2nes , 1973 ) Prinsenberg and Rattray, 1975 , Torgrimson 
52550 Hickey, 1979). 

In a brcader context, the oceanic response to hurricane 
Frederic was a geostrophic adjustment problen, which had 
both transient and stsady state components. After the tran- 
sient, inertic-qravity waves propagated away from the storm 
track, a gecstrophically kalanced ridge and current system 
remained under tke storm track. That is, the steady state 
currents were due to the balance between the horizontal 
pressure gradients and the Coriclis force (Geisler, 1970). 
Althcugh some adjustments to the mean flow are described, 
additional study is required to completely understand the 
interelaticnships between the forced mean flow, the forced 


wave field and the bottom topography. 
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AFPENDIX A 


A. CCMPLEX DENFODULAT ION 

The motivation fcr ccmplex demodulation is to isolate 
the carrier wave of a certain frequency, which must be kncwn 
a priori, and to form amplitudes and phases cf the modulated 
Signal. A linear filtering approach is used to form the 
instantaneous amplitudes and phases in contrasz to the har- 
monic analysis method. Essentially, the amplitude and phase 
time series represents a lccal harmonic analysis rather than 
discrete estimates averaged over continuous periods (Otnes 
and Enochscn, 1978). 

A band-limited time series is multiplied by the trigono- 
metric arguments of the carrier wave. The resultant time 
series is then lew-pass filtered to eliminate some high fre- 
quency noise that is generated. The cosine and sine coeffi- 
cients are ther combined tc form the amplitude and phase of 


the modulated wave. Consider the time series: 


x(i)=A(i) cos((2 7 fiat) + 0(1)) , (1) 


where fc is the carrier frequency, At is the sampling inter- 


val, @ is the phase , A is the amplitude of the wave, and i 
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represents a time index. Both the phase and amplitude vary 
as functicns cf time. Multiply equation (1) by the sine and 


cosine arguments of the carrier frequency f.e , 


x, à) =a (2) {cos (2 mf. iAt) + p(i))} sin (27 £. iAt) , (2) 
x, (3) 24 (i) {cos ((2 m 1At) + e (i0) cos (2 r £e “27007 (3) 
vhere 


Kel) c€0X(1) sin (2 r £f,1At) , and 


xod) = x (i) — (2526). 
Equation (2) can be revritten as 
X ç (i) 2 A (1)72 (SİR ((4 л22.24А2) + JS sin o (1) . (4) 
Similarly, equation (3) can be expanded and simplified as: 
70) = A(1)/2 {cos ((4 Tf. 1iAt) + en. COSTO) x5) 


Low pass filtering the abcve X4, ,x, series yields: 


Ys (i) =-a0)/2 [sin oc), (6) 


y, (i) A (1) /2 [cos $ ai 


where y. (1) and yc(i) are the low-pass filtered x,(i) and x 
(1). These data are filtered by convolving the input series 


with a set cf filter weights in the time domain (Bendat and 
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Ел=г=о1, 1971). Squaring and summing yo (7 and yu) 


yields 
asqı) sy + ye (2)) = Be 74 
er, 
Ali) = 2 a(i) (7) 
where A(i) is the amplitude as a function of tine. Simi- 


larly, the phase as a function of time is computed as: 


$(i) » Tanc! (-ус̧(2)/ 7, (2)) - (8) 


B. EOTENTIAL ENERGY 

The potential energy is computed to determine the rela- 
tionship of tke isotherz displacements to internal wave 
moticn. Given a temperature time series from a moored array 
and the ccrrespcnding vertical temperature qradien- df?/dz, 
the isotherm displacement ¿£ is computed in the following 
manner: 

7- —..-. ) > 
dr üz 

vhere T" is tle fluctuating part of the observed tempera- 


tures. The potential energy per unit volume is defined by: 
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p | 
PER ә ии (10) 


where N? is the Erunt-Vaisala frequency, ¢ is the isothern 
displacement fluctuations, and p, is the reference density. 


which is assumed to be unity. 


Ea ERAST SQUARES FIT 

Given a demodulated time series of the amplitude coeffi- 
cients and the vertical structure of the horizontal velocity 
eigenfunctions, the problem is tc compute the time-varying 
coefficients associated with each dynamical mode. The 
eme- varying atfplitudes U (Z; , t) code polist ezio, <... Za 


may be written as 


11 .2)= 5 С,(2)2, (21), 


Met 
where E,(z;) are the eigenfunctions for modes 1,2,3..., and 
Cy(t) are the ccefficients to be calculated in the analysis. 
Buc three tarcclinic  mscdal ccefficients are computed 
because current meters vere deployed at three depths in the 
water column. For example, the current at depth z,, is 


written as 
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and there are similar expressions at depth 2, and 23. The 
square of the error in the least squares fit at depth z, is, 
ee = to, - (C, Ej (Z4) + Cafa(Zi) + Cs Es (Z1 10} 2 (11a) 


and similarly at z and z , 


e2 = 10, - (COENA) + CoE, (zs) + Ca Fa (2,1) be , (11b) 


e$ 


iu, - (GEN. C CIE (z.) + CaE3 (2,1) 35 500 


The error is rinimized by solving the following set 


770 toe, /0C,+deg /0C,= 0 ; (12a) 
де2 /2С̧ *де2 /AC, +022 /0C,= 0 , (12b) 
SERIO oez oct oe? 70C,= 0 : (12c) 


Di fferentiating equation (11) with respect to C, and summing 


the expressions according tc equation (12a) yields 
U, E, (Z4) * U; E; (Z4) * U5E4(2,) = 
C, (E, (2,) E, (Z4) * E,(24) E, (22) + E, (23) E, (23)) 
+ C2(E, (2,) Ez (2,) + E, (22) EqQ(Z2) + E, (23) E2(23)) 


+ C4 (E, (24) Ey (2, ) + E,(Z2) E,(2,) + E,(Z3) E; (240); 


OI 


3 3 
Vara i m сс, (Zi) E,(zi) + 5 a) E 4 (Zi) 
s4 (al 


cal 
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3 | 
e BECOMES ZO (13) 


ira 


Now if equaticns (12b,c) are applied to equations (11), the 


corresponding equations for C, and C, become: 


3 3 
> F. (Zi) = 332 (zi) E,(Z,) + Co) B (zi) E, (Zi) 
xL "i "1 
f DINA E, (24) , (14) 
izi 
3 3 3 
Bau Enz) = C,) E, (20) ES (ei) + Co) 52171) E, (2;) 
i91, (x1 4=1 
3 
+ C4 E 4 (Zi) Ea (Zi) : (15) 
“L 


By inspection, let 


3 
Ann” En (Z:) En(Z;) TorBbum,uü-İl,Z.İI V (16) 


izi 
ə 


Bu =) Ui Eulzi) (17) 


i1. 


Therefore, equations (13-15) Gan be simplified using tha 


above definiticns inte a simple matrix =quation of the form 
(18) 


Equation (18) isa simple symmetric matrix equation and it 


was sclved using the IMSL routine LEQT2F on the IBM 3033. 
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AEPENDIX IE 


A. DEEP REGICN 

The ocean currents at CMA2,3 during the quiescent period 
clearly showed topographical effects associated with the 
DeSotc Canyon. The behavicr of the surface currents is best 
described by a EVD (Fig. 29). The strong u component at 
CMA2 produced a total displacement of 600 km over 39 days 
for a 15.4 kmyd mean current. Inertial waves were super- 
posed on this trean flow, particularly during the period from 
555 2 to 239. The mean flow increased during this period. 
At CMA3, the mean flow decreased in the southerly direction 
5 0:1..5 to 9.5 km/d (Fig. 29b). The different orientation 
in the mean flow across the entire shelf region is presun- 
ably due to the influence cf bottom topography. 

The trajectory of the near-bottom horizontal currents at 
CMA2 (Fig. 30a) was in the same direction as in the mixed 
Mayer until JD 222. After the reversal, “¿he mean flow was 
towards the west at about 3.3 knmyd. The bottom PVD at CMA3 
(Fig. 305) showed ttat iritially the mean flow was towards 


the north at 2.0 km/d over the first 6 IP. The £low then 
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Figure 29. Surface PVD fer a) CMA2 5: b) CMA3. ə i 
represents a time span 48 start 5. 
n for CMA2 (19 m) and . JD 213 for "V 
m). 


0. 


Changed direction and accelerated tc about 2.8 km/d in the 
same directicn as in the surface layer. 

Horizontal kinetic energy spectra for this period are 
summarized in Table XI for both inertial and semi-diurnal 
frequency tands. The vertical variation in the inertial HKE 


levels at CMAZ indicated a larger amount of energy in tha 
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27 0 30. Bettcm PVD for a) CMA2 and b) CMA3. The = 
represents a time span of 48 h starting a) 
“I for CMA2 (324 m) and b) JD 213 for сиз" (457 
m). 
thermocline than at the surface or near-botton. However, 
the difference ketween the surface and bottom energy levels 
vas abou- cne crder of magnitude. At CMA3, the inertial HKE 
estimates decreased steadily with depth and there was more 


than an crder of magnitude difference between the surface 


and bottom layers. The HKE levels were slighty greater at 
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CMA3 in tke surface layer and thermocline, but they were 


slighty less near the bottcn than at CMA2. 


TABLE XL 
Normalized HKE Spectral Estimates at 957 Confidence 


Normalized 


Meter Nos Of Data . HKE Spectra 
Depth Points Bandwidth M2 
(m) cph (cm/s) t27cph 
a) CMAZ 
19 5780 0.001 MATOS ee 
179 4590 0.0013 5727710777 0XXx 001 
324 5780 05007 2277277” “0 О 10 ә 
b) CMA3 
21 4680 05000: 6557507770 0” 
251 4680 0.0073 4.5x10+3 7.0x10+2 
437 4680 0.0013 TES IOS "476X1075 
457 4680 0.0013 2. 0x 1082 1.7x10t+t1 


* instrumentaticn prcblens with time clock 
M2 is the semi-diurnal tide frequency band 
f is the inertial/diurnal tide frequency band 

The semi-diurnal HKE estimates also varied spatially. 
The estimates cf HKE in tte thermocline and near-bottom lay- 
ers at CMA2 were not significant; however, “he HKE estimates 
at CMA3 were significant at all depths. The HKE estimates 
at 437 m even exceeded both the surface layer and thermo- 


cline levels ty about an czder of magnitude. In comparison 


with the near-bottom estimates at 457 m, the kinetic energy 
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levels were mcre than two orders of magnitude greater. 
These proncunced changes in “the semi-diurnal energy levels 
with depth indicate the presence of internal tide motion. 
In most of tte kinetic energy spectra examined, longer 
pericd motions were insignificant, because their energies 
were smeared acrcss adjacent bands except in the bottom lay- 
ers. Subinertial variability at periods of about 3 IP was 
quite apparent in the bottcm reccrds at both CMA2 and CMA3. 
These moticns were generally more energetic than the semi- 
diurnal tides. 

Rctary spectra are analyzed fer the surface and bottom 
Current records and are given in Table XII. Most of the HKE 
2 otility is ccnsistent vith  pclarized, CW rotating iner- 
214) motions. The CW energy level exceeds the CCW level by 
about an crder of magnitude. The ellipse orientations for 
these waves retates CCW (kacks) with depth. These ellipses 
are stable except for the near-bottom inertial period motion 
at GMA2. The direction of the semi-major axis of the 
ellipse varies over time, which indicates that the inertial 
waves are isotropic. The semi-diurnal tidal motion is 
polarized in tke CW directicn. Generally, the semi-diurnal 


tidal moticns are coherent, although the phase differences 
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between the hcrizcntal velccity components vary between the 


arrays. 
TABLE XII 
Ncrmalized Rotary Spectrum Analysis 
Normalized 
Rotary 
Spectra Ellipse | 
ER Depth CH CCW Stap: Dir. Rotary 
Bang (2) (cmZs) +2/cph (Deg-) Coeff. 
a) CHAZ 
I 19 1.6x10+3 1.9x10+2 0.80 103 5570 
S 19 1.6x190t2 3.6x10*! 0.39 101 +0,87 
T 324 3.0x10* 2 1.4x10+1 0.05 90 +0.37 
5 324 1.0219 9 11x 10807 0.60 169 059 
b) CMA3 
m 21 1. Zx10** 1.7x 10+ 1 0.40 157 + Ol a 
S 21 4.0x10+2 3.0x 10+ 1 0.61 DD +0.63 
T 457 2-10“ 1.6x10*1 Qe. | 138 0.986 
S 457 3.022081 3.9x10*9 0.20 93 EOS 


I=inertial frequency ba nd 
S=Semi-diurnal frequency tand 


Ocean kottcn temperatures ranged between 8.0 and 12. 


22-X-0:2 and 8.0 to 9.90 0Ç at CHA3 (Fig. 31), but there 
different variations and trends. At CMA2, (not shown) 
temperature decreased gradually to about 8.3 ОС near JD 
and then increased abruptly withcut any periodicities. 


tom temperatures at CMA3 were much lower and “here 


SOC 
vera 

the 
2407 
Bet = 


Was 


near-inertial motion superfosed on this early cooling trend. 
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The minimum temperature occurred on JD 232, which was about 
2 IP tefore tte rinimum otserved at CMA2. These minima in 
the temperature, in additicn to the cooling trends, are man- 
ifestations of the advection of cold water from +he DeSotc 


Carycn by tte mean 21сҹ. 


B. SHALLOW REGICN 

In the thermocline (49 m), the inertial oscillaticns 
were not nearly as energetic and did not rotate CW as in «he 
mixed layer (Fig. 32): The non-stationary trends ir the 
east-west currents were similar to the mixed layer observa- 
*icns, but tkere was an initial mean flow towards the deep 
ocean. Superposed on the mean flow were lower frequency 
ouiesttatiens with periods cf the crder of about 10 to 14 IP. 
The total displacement over the forty days vas about 300 km 
5.:ch correspondad to a mean current of 7.5 kn/d. 

Inertial cscillations at 64 m were not as energetic as 
at 49 m, but there was a similar non-stationary trerd in the 
time series. This  *rend was associated with a lower fre- 
quency waves having a pericd of 10 to 14 IP. These oscilla- 
tions were superposed on a northward mean flow which was 
slighly less than 7.0 km/d. Records from the current meter 


nearest to the bottom (92 m), were out of synchronization 
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207üre 32. Mid-depth PVL at CMA1. The * represents a time 
span of 48 h starting on JD 207 at 49 m. 
due to internal clock problems. For example, large 
increases in the velocity cccurred on JD 252 as opposed to 
256. Hence, tte data were not included. 
The HKE spectra for the 21, 49, and 64 m depths are 
shown in Table XIII for bcth the inertial and semi-diurnal 


frequency Lands. The inertial/diurnal tidal kinetic energy 
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in the surface layer was two orders of magnitude larger than 
in *he sutsurface layers. The semi-diurnal tidal motion 
decreased from the surface to 49 m and then increased with 
depth such that the HKE near the bottom almost equaled that 


near the surface. 


TABLE KILL 
Ncrmalized HKE Spectral Estimates at 95% Confidence 


Normalized 


Meter No. of Data Frequency HKE Spectra 
Depth Points Re sclution f M2 
(2) cph (cn/s)*?/cph 
CHAI 
21 5825 0.001 1.4x10+* 1.5x10+2 
49 625 0001 1.4x10+3 4.3x10+1 
64 5825 005001 5.0x10+2 8.0x10+1 
92% 5825 0.001 1.5x10+2 1.0x10+2 


* instrumentaticn prcblens with time clock 
M2 is the segi-diurnal tide frequency band 
f is the inertial/diurnal tide frequency band 

Temperature data at kE cth 21 and 49 m are not available 
because of the thermistor problem. The temperatures ranged 
maomeio.09°C to the thermistor cutoff of 21.5 OC at 49 m. 
Temperature variations were about 0.59C over an IP, although 
most of the variability was associated with some non-sta- 


tionary trends and lower frequency motions similar to those 


observed in tbe currents. Just prior to the storm, tempera- 
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tures decreased to a relative minimum as cold water was 
advected Ey the mean flow. 

In summary, based on the observations during the quies- 
cent pericd, the mean flew strcngly depends on the bottom 
topography. Inertial band motion is evident at all depths 
during the cuiescent pericd. Scone of chis notion is due tc 
freely prcpagating inertial waves while the remainder is 
Bes of the fcrced diurnal tides. Internal tides, as well 
as inertial oscillaticns, are part of the freely propagating 
internal wave continuum below the surface layer. The effect 
of the topcgraphy on the mean flow and the wave motions fur- 
ther complicates the dynamics of the region and the resul- 
25:55 circulation. Longer period waves (2.5 to 3 day 
pericd) set ur north-south oscillations that advect cooler 
water from the DeSoto Canycn onto the shelf region. Super- 
inertial frequency motions are influenced by the diurnal and 


semi-diurnal tides and contribute to the variability. 
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